THIS  REPORT  HAS  BEEN  DELIMITED 
AND  CLEARED  FOR  PUBLIC  RELEASE 

UNDER  DOD  DIRECTIVE  5200,20  AND  \ 

NO  RESTRICTIONS  ARE  IMPOSED  UPON 
ITS  USE  AND  DISCLOSURE, 

DISTRIBUTION  STATEMENT  A 

APPROVED  FOR  PUBLIC  RELEASE; 


DISTRIBUTION  UNLIMITED, 


DEFENSE  DOCUMENTATION  CENTER 

FOR 

SCIENTIFIC  AND  TECHNICAL  INFORMATION 

CAMERON  STATION.  ALEXANDRIA.  VIRGINIA 


UNCLASSIFIED 


THIS  DOCUMENT  IS  BEST 
QUALITY  AVAILABLE.  THE  COPY 
FURNISHED  TO  DTIC  CONTAINED 
A  SIGNIFICANT  NUMBER  OF 
PAGES  WHICH  DO  NOT 
REPRODUCE  LEGIBLY. 


NOTICJ:  Vhen  govertaent  or  other  drawings,  speci¬ 
fications  or  other  date  axe  used  for  any  purpose 
other  than  In  connection  vith  a  definitely  related 
government  procurement  operation,  the  U.  S. 
Qovsrnaent  thereby  incurs  no  responsibility,  nor  any 
obligation  whatsoever;  and  the  fact  that  the  Qove re¬ 
sent  aay  have  formulated,  furnished,  or  in  any  way 
supplied  the  said  drawings,  specifications,  or  other 
data  is  not  to  be  regarded  by  implication  or  other¬ 
wise  ae  in  any  Banner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights 
or  permission  to  Manufacture,  use  or  sell  any 
patented  invention  that  aay  in  any  way  be  related 
thereto. 


6  H 


cc 

•  o> 

© 


* 

02 

>, 

Tt* 

7J 

*  -s 

» 

p  2 

o 

Li— 

CT- 

C' 

. 

'-■C 

c 3 

<C 

NAVWEPS  REPORT  8442 
NOTS  TP  3399 
COPY  59 


A  LIFTING-SURFACE  PROPELLER  DESIGN  METJ1D0 
F93  HIGH-SPEED  COMPUTERS 

by 

D.  M.  Nelson 

Underwater  Ordnance  Department 

Released  to  ASTIA  for  further  dissemination  v;ith 
out  limitations  beyond  thoso  iupoued  by  security 
regulations. 


ABSTRACT.  A  propeller  design  method  that  treats  the 
blades  as  lifting  surfaces  has  been  developed  and  pro¬ 
grammed  on  an  IBM  7090  computer.  The  singularity 
distributions  representing  the  bound  circulation,  shed 
vorticity,  and  blade  thickness  are  treated  as  continuous 
sheets,  making  it  necessary  to  determine  the  Cauchy 
principal  value  of  the  improper  integrals  expressing 
the  induced  velocities  on  the  blades.  The  camber  lines 
and  pitch  angles  are  determined  from  the  variation  of 
the  normal  component  of  induced  velocity  across  the 
blade.  This  method  is  applicable  to  moderately  loaded, 
non-optimum,  wake-adapted  propellers  with  or  without 
skew. 
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FOREWORD 

This  report  presents  the  development  of  a  propeller  design 
method  that  treats  the  blades  as  lifting  surfaces.  The  method 
has  been  programmed  on  an  IBM  7090  computer  and  is  available 
for  design  calculations. 

The  work  was  undertaken  to  make  available  at  the  U.  S.  Naval 
Ordnance  Test  Station  the  most  up-to-date  technique  for  the  de¬ 
sign  of  marine  propellers.  It  was  done  from  September  1961  to 
October  1963  under  Bureau  of  Naval  Weapons  Task  Assignment 
RUTO-3E-000/216-1/R009-01-03,  Problem  Assignment  428. 

The  considered  opinions  of  the  Propulsion  Division  are 
represented  in  this  report. 
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NOMENCLATURE 


Chord 


Chordwise  coordinate  measured  from  leading  edge 
(Fig.  3a) 

Infinitesimal  area  of,  helical  sheet 
Infinitesimal  length  of  vortex  line  element 
Propeller  diameter 
V  C 

- rR 


fi»f?» 


A 


Functions  of  x  (Eq.  10) 

Number  of  blades 

Nondimensional  bound  circulation 

Nondimcnsional  camber  offset  (Fig.  13) 

Ratio  of  the  strength  of  a  line  element  of  'vee  vorticity 
at  any  chordwise  station  to  that  at  the  traii.‘ ng  edge 

Distance  from  the  stacking  line  to  the  leading  edge  of 
blade  section  (Fig.  3a) 

Radial  coordinate 

xatio  of  percent  thickness  of  b.’ade  section  at  any  span- 
wise  station  to  percent  thickness  at  reference  station 

Propeller  radius 

Distance  from  a  point  on  one  of  vhe  blades  to  point  p 
Local  axial  inflow  velocity  to  prope.Ue* 

Forward  speed  of  propeller  relative  to  undisturbed  fluid 
Relative  velocity  of  blade  section  and  flui'-  (Fig.  12) 
Induced  velocity 

Normal  component  of  induced  velocity 

Contribution  to  from  the  singularity  region 

Nondimensional  radial  coordinate  (Eq.  7) 

Nondimensional  radial  coordinate  measured  from  center 
of  chordwise  strip  or  center  of  singularity  .-egion 
(Fig. 7,  9,  11) 
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x,  y,  z 

x',  y z* 
X',  Y',  z* 
y 


Nondimensional  rectangular  coordinate  system  at  pro¬ 
peller  axis  (Eq,  7) 

Rectangular  coordinate  system  at  propeller  axis  (Fig.  1) 

Rectangular  coordinate  system  on  helical  sheet.  (Fig.  2) 

Nondimensional  chordwise  coordinate  (Eq.  4) 

Nondimensional  chordwise  coordinate  measured  from 
center  of  spanwise  strip  or  center  of  singularity  region 
(Fig.  4,  7) 

Nondimensional  half- thickness  of  blade  (Fig.  12) 


a  Angle  defining  points  on  helical  sheets  (Fig.  1) 

a ^  Angle  defining  points  of  helical  sheets  measured  from 

center  of  singularity  region 

P  Pitch  angle  of  helical  sheets  (Fig.  1) 

y  Variation  of  pitch  angle  from  lifting-line  value 

T  Bound  circulation 

T  Strength  of  vortex  line  element 

q  Source  strength 

X  x  tan  p 

p  Surface  source  density 
<}>  Source  potential 

ijj  Angle  describing  stacking-line  locations  (Fig.  1) 
m  Angular  velocity  of  propeller 


SUBSCRIPTS 

B  Due  to  bound  circulation 

F  Due  to  free  vorticity 

p  Corresponding  to  the  singularity  point  p 
T  Due  to  blade  thickness 
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INTRODUCTION 

The  existing  marine  propeller  design  methods  that  can  be  utilized 
without  a  high-3peed  computer  are  normally  built  up  from  optimum 
propeller  lifting-line  theory,  two-dimensional  slender-airfoil  theory, 
and  correction  factors  applied  to  the  two-dimensional  results  to  account 
for  three-dimensional  effects.  Of  these  methods  the  most  comprehen¬ 
sive  is  probably  that  of  Eckhardt  and  Morgan  (Ref.  1). 

Since  these  methods  are  applied  to  non-optimum  propellers  of  es¬ 
sentially  arbitrary  blade  shape  and  even  to  wake-adapted  propellers, 
close  examination  of  their  underlying  aosumptions  raises  some  ques¬ 
tion  about  their  ability  to  produce  satisfactory  propeller  designs.  For 
example,  the  induced  velocities  at  the  lifting  line  predicted  by  optimum 
propeller  theory  are  positive  at  all  spanwise  stations  on  the  blade. 
However,  as  Lerbs  (Ref.  2)  points  out,  and  as  is  substantiated  by  the 
present  work,  non-optimum  propellers  may  and  normally  do  have 
rather  large  negative  induced  velocities  near  the  hub  and  the  tip.  Also, 
the  camber  correction  factors  that  usually  have  been  obtained  from  the 
work  of  Ludwieg  and  Ginzel  (Ref.  3  and  4)  are  applied  to  both  free- 
running  and  wake-adapted  propellers  of  arbitrary  blade  shape,  circu¬ 
lation  distribution,  and  advance  ratio.  Since  the  Ludwieg  and  Ginzel 
calculations  were  made  for  only  a  few  cases  of  free- running  propellers, 
rather  crude  interpolations  on  extrapolations  are  necessary.  Further¬ 
more,  the  pitch  correction  factor  proposed  by  Lerbs  (Ref.  5)  and  used 
in  the  Eckhardt  and  Morgan  method  is  derived  on  the  assumption  that 
the  propeller  is  developing  lift  by  angle  of  attack  and  is  then  applied  to 
propellers  developing  lift  by  camber.  Inconsistencies  such  as  these 
may  be  the  cause  of  the  inadequacies  of  these  methods,  because  ex¬ 
perience  has  shown  that  they  do  not  always  lead  to  satisfactory  pro¬ 
peller  designs,  particularly  when  cavitation  resistance  is  concerned. 

Lack  of  confidence  in  such  approximate  methods  has  prompted  a 
considerable  amount  of  work  in  the  past  few  years  on  lifting -surface 
theory  and  design  methods  for  marine  propellers.  Work  in  the  United 
States  has  been  done  by  Pien  (Ref.  6)  and  Kerwin  (Ref.  7),  in  the  Nether¬ 
lands  by  Sparenberg  (Ref.  8)  and  by  van  Manen  and  Bakker  (Ref,  9),  in 
England  by  Cox  (Ref.  10},  and  in  Japan  by  Nishiyama  and  Nakajima 
(Ref.  11).  In  all  these  studies,  however,  no  extensive  calculations 
were  made  that  would  permit  the  design  of  a  non-optimum,  wake- 
adapted  propeller  with  a  specified  blade,  sjiape  and  circulation  distri¬ 
bution.  In  fact,  the  number  of  parameters  that  would  have  to  be  varied 
(wake  fraction,  blade  shape,  circulation  distribution,  number  of  blades, 
and  advance  ratio)  and  the  long  computing  times  needed  for  such  calcu- 
tions  make  such  an  investigation  impractical.  Hence,  to  design 
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propellers  treating  them  as  lifting  surfaces,  the  facilities  must  be  at 
hand  to  carry  out  the  design  calculations.  To  this  end,  the  work  de- 
scilbc«Pirrih^^.r«po'rt*wanPurTdwtakwnV’'  —  ' 

Propellers  for  torpedoes  are  of  primary  interest  at  the  Naval 
Ordnance  Test  Station.  Since  these  propellers  operate  in  a  wake,  and 
since  good  cavitation  resistance  dictates  circulation  distributions  that 
are  not  optimum,  a  lifting- surface  design  method  that  is  applicable  to 
non-optimum,  wake-adapted  propellers  has  been  developed.  This  de¬ 
velopment  goes  one  step  beyond  any  of  the  works  referenced  above  in 
that  the  effect  of  the  blade  thickness  is  accounted  for. 


REPRESENTATION  OF  PROPELLER  SYSTEM 
BY  SINGULARITY  DISTRIBUTIONS 

Each  propeller  blade  is  represented  by  a  continuous  sheet  of  vor- 
ticity  cre&img  the  bound  circulation  and  by  a  continuous  sheet  of 
sources  and  sinks  representing  the  blade  thickness.  Since  elements 
of  vorticity  cannot  end  in  space,  spanwise  variations  in  the  bound  cir¬ 
culation  give  rise  to  elements  of  free  or  shed  vorticity  in  the  fluid  that 
are  left  in  helical  paths  along  and  behind  the  propeller.  The  shed  vor¬ 
ticity  for  each  blade  is  represented  by  a  continuous  helical  sheet  of 
vorticity  starting  at  the  leading  edge  of  the  blade  and  extending  back  to 
infinity.  In  this  development,  the  hub  boundary  condition  has  been  ig¬ 
nored  and  there  are  therefore  no  singularities  representing  the  hub. 


SPATIAL  LOCATION  OF  SINGULARITIES 
REPRESENTING'  PROPELLER  SYSTEM 

Before  the  calculation  of  the  velocities  induced  by  the  propeller  can 
begin,  the  singularities  representing  the  propeller  system  must  be 
properly  located  in  space.  Lerbs*  induction  factor  method  (Ref.  2)  for 
the  lifting-line  solution  of  the  moderately  loaded,  non-optimum,  wake- 
adapted  propeller  offers  the  best  approximation  to  the  spanwise  pitch 
distribution  of  the  helical  sheets  upon  which  the  singularities  represent¬ 
ing  the  propeller  should  be  located.  In  fact,  for  the  free-running,  un¬ 
skewed,  optimum  or  lightly  loaded  non-optimum  propeller  with  negli¬ 
gible  blade  thickness,  symmetrical  blade  shape,  and  symmetrical 
chordwise  loading,  the  spanwise  pitch  distribution  given  by  Lerb3' 
lifting-line  solution  is  identical  to  the  spanwise  pitch  distribution  of  the 
blade  chordlines  given  by  the  lifting -surface  solution.  For  other  condi¬ 
tions,  the  pitch  of  the  blade  chordlines  will  usually  differ  only  slightly 
from  the  pitch  of  the  helical  sheets  given  by  Lerbs*  lifting-line  solution. 

Since  the  object  of  the  liftirtg- surface  solution  is  to  determine  the 
camber  lines,  it  is  not  possible  to  distribute  the  vorticity  representing 
the  bound  circulation  on  the  camber  lines.  Instead,  the  vorticity  is 
distributed  on  the  helical  sheets  from  Lerbs*  solution.  This  approxi¬ 
mation  is  similar  to  the  one  made  in  two-dimensional  slender- aix*  foil 
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theory  in  which  the  circulation  is  distributed  along  the  chordline  when 
the  camber  lines  are  computed.  Since. o.uly  a.TO2ll..jxi3QUjaSi-.o£  csmbex.  *• 
are  needed  in  marine  propellers  to  generate  the  desired  lift,  this  ap¬ 
proximation  can  lead  to  no  serious  error  in  computing  the  camber  lines. 
The  sources  and  sinks  representing  the  blade  are  also  distributed  on 
the  helical  sheets  when  the  effects  of  blade  thickness  are  computed. 

Since -the  free  or  shed  vorticity  in  the  fluid  is  left  behind  in  helical 
paths  whose  pitch  is  determined  by  the  relative  flow  of  the  propeller 
and  fluid,  the  elements  of  free  vorticity  fall  naturally  on  helical  sheets 
that  are  very  close  to  or  identical  to  those  given  by  Lerbs'  lifting-line 
solution.  Hence,  when  the  effect  of  the  free  vorticity  is  computed,  it 
is  distributed  on  the  helical  sheets  from  Lerbs'  solution. 


COORDINATE  SYSTEM?  USED  IN  ANALYSIS 

Figure  1  shows  the  coordinates  used  to  describe  the  helical  system 
upon  which  the  singularities  are  distributed.  The  angles,  ijjm,  which 
are  given  by 


m  -  1 

+m  ~ - 2tt  (m  =  1,  2,  . .  .  ,  g) 

g 

where  g  is  the  number  of  blades,  describe  the  lines  where  the  propeller 
blades  intersect  the  plane  X'  =  0.  These  lines  are  referred  to  as  the 
stacking  lines  and  are  chosen  coincident  with  the  lifting  lines  of  Lerbs' 
solution.  For  unskewed  propellers  they  will  also  be  chosen  coincident 
with  the  line  passing  through  the  centers  of  lift  of  the  blade  sections. 


Z' 
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One  line  element  at  a  constant  radius  of  the  helical  sheet  streaming 
back  from  one  of  the  propeller  stacking-line  locations  is  shown.  This 
element  makes  the  pitch  angle  p  with  respect  to  the  plane  X'  =  0.  The 
coordinates  of  a  point  on  the  helical  sheets  in  the  rectangular  X',  Y\  Z1 
coordinate  system  can  then  be  written 

X'  =  ar  tan  p 

Y'  =  -r  ain(4<m  +  a) 

Z'  =  r  cos(4>m  +  a)  (1) 

Another  coordinate  system  needed  in  the  analysis,  the  rectangular 
X1,  Y',  Z'  system,  is  shown  in  Fig.  2.  In  this  system  the  Z'  axis  inter¬ 
sects  the  X'  axis  perpendicularly,  the  X'  axis  is  tangent  at  point  p  to  a 
constant  radius  line  element  of  the  helical  sheet  streaming  back  from 


the  stacking-line  location  ipl  =  0,  and  the  Y'  axis  is  normal  to  this  line 
element  at  point  p.  A  vector  having  components  Ax,  Ay,  _Az  in  the 
X',  Y'»  Z'  coordinate  system  will  have  components  in  the  X',  Y',  Z' 
coordinate  system  given  by 

Ax  =  Ax  sinPp  -  Ay  coscip  cosPp  -  A^  sinap  cos  Pp 

Ay  =  Ax  cosPp  +  Ay  cosap  sinPp  -f  Ajr  sinap  sinpp 

Az  =  -Ay  sinap  +  Ajj  cosap  (2) 
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APPROACH  TAKEN  IN  MATHEMATICS  OF  SOLUTION 

When  the  velocities  induced  oy  the  singularities  representing  the 
propeller  system  are  computed,  double  integrals  must  be  solved.  The 
integrands  are  such  that  integration  cannot  be  performed  analytically 
with  respect  to  either  of  the  variables.  In  order  to  bypass  the  diffi¬ 
culty  of  performing  numerical  double  integration,  the  region  to  be  in¬ 
tegrated  is  divided  into  strips  sufficiently  narrow  that  certain  approxi¬ 
mations  hold  accurately  across  them.  These  approximations,  together 
with  some  restrictions  put  on  the  form  of  the  circulation  distribution, 
allow  integration  to  be  performed  across  the  strips  analytically.  Thus, 
instead  of  a  numerical  double  integration  there  is  an  analytical  integra¬ 
tion  across  the  strip,  a  numerical  integration  along  the  strip,  and  a 
summation  over  the  strips.  This  process  is  applied  everywhere  except 
in  a  small  region  around  the  point  where  the  velocity  is  desired,  since 
a  singularity  occurs  in  the  integrand  at  that  point,  making  numerical 
integration  through  the  point  impossible.  Integration  is  carried  out 
analytically  over  this  small  singularity  region  to  obtain  the  Cauchy 
principal  value.  To  make  an  analytical  integration  possible  over  this 
region,  several  approximations  were  necessary  which  in  essence 
amount  to  (1)  assuming  a  linear  variation  of  most  of  the  variables  over 
the  region,  plus  (2)  discarding  third  and  higher-order  terms  in  the 
small  nondimensional  coordinates  describing  the  distance  of  a  point 
from  the  singularity  point.  The  check  solutions  run  to  date  showed 
that  the  region  may  be  chosen  sufficiently  small  that  these  approxima¬ 
tions  yield  a  good  value  of  the  contribution  of  the  region  to  the  velocity, 
without  being  so  small  that  the  numerical  integration  in  the  area  im¬ 
mediately  surrounding  the  region  cannot  be  carried  out  accurately. 

An  important  feature  of  this  approach  is  that  the  point  where  the 
velocity  is  desired,  i.  e.  ,  the  singularity  point,  is  always  chosen  at 
the  center  of  the  strips  and  at  the  centroid  of  the  small  region  sur¬ 
rounding  tlic  singularity.  Also,  all  approximations  are  expanded  about 
the  center  of  the  strips  or  the  centroid  of  the  singularity  region.  By 
doing  this,  the  accuracy  of  the  method  is  greatly  enhanced  and  comput¬ 
ing  times  are  reduced  since  the  number  of  points  needed  to  obtain  an 
accurate  numerical  integration  is  smaller. 

Two  restrictions,  mentioned  above,  were  placed  on  the  circulation 
distribution  to  facilitate  integration  across  the  strips:  (1)  the  chord- 
wise  distribution  of  circulation  must  be  made  up  of  straight  line  seg¬ 
ments,  and  (2)  the  spanwise  distribution  of  circulation  must  lend  itself 
to  an  accurate  piecewise  approximation  by  parabolic  sections.  Since 
fairly  optimum  chordwisc  distributions  of  circulation  from  the  stand¬ 
point  of  cavitation  resistance  can  be  made  up  of  straight  line  segments, 
and  since  most  practical  spanwise  circulation  distributions  do  not  have 
extremely  large  third  derivatives,  these  restrictions  offer  no  serious 
restriction  on  the  generality  of  the  method. 
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DESIRED  COMPONENT  OF  INDUCED  VELOCITY 

Because  the  component  of  induced  velocity  parallel  to  the  blade- 
section  chordline  will  normally  be  very  small  compared  to  the  relative 
velocity  of  the  blade  section  and  the  still  fluid,  it  may  be  neglected  in 
the  boundary  condition  when  the  camber  lines  are  being  determined. 
Only  the  component  of  induced  velocity  normal  to  the  blade- section 
chordline  is  needed  to  determine  the  camber-line  boundary  condition. 
This  approximation  is  identical  to  that  used  in  slender-airfoil  theory. 
Assuming  that  all  induced  velocities  are  to  be  computed  on  the  blade 
i|jj  =  0,  the  normal  component  of  induced  velocity,  denoted  by  Wy,  will 
be  in  the  Y'  direction  shown  in  Fig.  2.  There  are  three  contributions  to 
Wy:  the  one  resulting  from  the  bound  circulation  (Wy)B»  the  one  re¬ 
sulting  from  the  free  vorticity,  (Wy)p,  and  the  one  resulting  from  the 
blade  thickness  (Wy)y.  Hence, 

Wy  =  {Wy)B  +(Wy)F  +  (Wy)T  (3) 


CALCULATION  OF  NORMAL  COMPONENT 
OF  INDUCED  VELOCITY 

VELOCITY  DUE  TO  BOUND  CIRCULATION  (Wy)B 

The  position  of  points  on  a  blade  section  are  shown  in  Fig.  3a,  where 
C  is  the  chord,  d  is  the  distance  from  the  leading  edge,  and  L  is  the 


FIG.  3a. 

distance  from  the  stacking  line  to  the  leading  edge.  A  nondimensional 
chordwise  coordinate  is  defined  as 


d 

y  =  -  (4) 

c 

The  location  of  each  vortex  line  element  of  the  bound  circulation  sheet 
is  given  by  the  condition  y  =  Constant,  as  illustrated  in  Fig.  3b  for  an 
unskewed  propeller  and  in  Fig.  3c  for  a  skewed  propeller.  This  choice 
of  the  location  of  an  element  of  the  bound  circulation  sheet  is  the  most 
convenient  one,  since  elements  of  free  or  shed  vorticity  arise  only 
from  a  change  in  the  spanwise  distribution  of  circulation.  For  the 
other  methods  of  locating  the  elements  of  the  bound  circulation  sheet, 
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(Blade  shown  flat  here  and  in  follow¬ 
ing  .sketches  to  simplify  drawing) 


Stacking  line 

Vortex  line 
element  of 
bound  cir¬ 
culation 
sheet 


y  =  Constant 

Unskewed  Propeller 
FIG.  3b. 


Stacking  line 

Vortex  line 
element  of 
bound  cir¬ 
culation 
sheet 


y  =  Constant 

Skewed  Propeller 
FIG.  3c. 


such  as  that  shown  in  Fig,  3d,  elements  of  shed  vorticity  arise  not  only 
from  a  change  in  the  spanwise  distribution  of  circulation  but  also  as  a 
result  of  blade  taper.  The  added  complication  of  having  to  consider 
both  these  effects  in  the  calculations  involving  the  free  vorticity  make 
such  a  method  impractical. 


Vortex  line 
elements  of 
bound  cir 
culation 
sheet 


Impractical  Representation 
of  Bound  Circulation 


FIG.  3d. 


Referring  to  Fig.  1  and  measuring  the  distances  d  and  L  along  the 
constant  radius  helical  line  with  the  stacking  line  in  the  plane  X'  =  0, 
the  following  relation  is  obtained  between  a  and  y  and  r. 

C  cos  p  L  cos  p 

a  =  y - 

r  r 
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Letting  the  subscript  p  denote  the  point  p  where  the  normal  compo¬ 
nent  of  velocity  is  desired,  and  remembering  th.ai,,aU. , valacltlaa-  are,  to„^.. 
be  computed  on  the  blade  =  0,  the  coordinatesrirr  the-  X',  Yr,  Z"r  sys¬ 
tem  of  the  point  p  are,  from  Eq.  5, 


xp  =  Yp(c  8inP)r=rp  “  (L  8inP)r=rp 

Yp  =  _rp  sinQp 

Zp  =  rp  cos Qp  (6) 

where 


(C  cos  p 


ap  =  yP 


f  L  cos  p 


r=rr 


'r=rT 


The  X',  Y',  Z*  components  of  the  distance,  S,  from  an  element  of 
vorticity  on  one  of  the  blades  to  the  point  p  are  obtained  from  Eq.  5 
and  6. 


Sx  =  Xp  -  yC  sinp  +  L  sinp 

(C  cos  p  L  cos  p 

+m  +  y - 

r  r 

1C  cos  p  L  cos  p 

+m  +  y - 

r  r 

Nondimensionalizing  with  the  propeller  radius,  R,  and  diameter,  D, 
and  introducing  the  nondimensional  coordinates, 

r  X'  Y1  Z' 


x  = 


R 


(?) 


the  components  of  di  and  S  may  be  written 


dix 

-  =  yf5  -  f6 

R.dx 


dfy 

- =  -cos  (ijjm  +  yf3  -  f4) [y(f 7  -  (3)  -  (f 8  “  4)1  ~  sin  (^m  +  y*3  "  4) 

Rdx 


-  =  -sin  {1 {/m  +  yf3  -  4)[y(f7  “  4)  “  ({8  ~  4)1  +  cos  (4>m  +  y4  “  4) 

Rdx 
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SX 

- =  Xp  -  yfi  +  f2 

R  p 

SY 

- -  Y p  +  x  aln  (i| <m  +  yf3  -  £4) 

R  r 

—  =  Zp  -  X  cos  (4*m  +  yf3  -  f4)  (8) 

R 

where 

Xp  =  Ypflp  “  f2p 
Yp  =  -xp  ainap 

Zp  =  Xp  COSQp 

Qp  =  Ypf3p  ‘  f4p  <9) 

and  where 


C 


{ l  =  Z  —  sin  P, 

D 

flp  =  (fl)x=xp 

L 

£2  =  2—  sin  p, 

D 

f2p  =  <f2)x=xp 

2{C/D)  cos  p 
f3  =  * 

X 

f3P  =  tf3>x=xp 

2(L/D)  cos  p 

f4  = - - 

X 

f4p  =  (f4)x=Xp 

d[2<C/D)  sin  p] 

f5  = - , 

dx 

^5p  =  ^5)Xr;Xp 

d[2(L/D)  sin  p] 

f6  = - - , 

dx  , 

f6p  =  <f6)x=Xp 

d[2(C/D)  cosp] 

f7  = 

dx 

f7p  =  (f7)x=Xp 

d[2(L/D)  cosp] 

f8  =  * 
dx 

f8p  =  (f8)x=Xp 

1C 


sw* 
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Breaking  the  blade  into  spanwise  strips  bounded  by  lines  y  -  Constant, 
as  illustrated  in  Fig.  4,  there  may  be  written  for  each  strip 

r**rn  +  7a’  ■  - 

The  subscript  n  takes  on  the  values  1,  2,  3,  ...  ,  s,  where  s  is  the 
number  of  strips.  Substituting  Eq.  11  into  Eq.  8,  using  the  trigonometric 


yi 

n 

Y3 


FIG.  4. 


\ 

\ 


\ 


I 


identities  for  the  sin  and  cos  of  the  sum  of  two  angles,  and  re¬ 
stricting  the  maximum  half-width  of  the  strips,  (y^)max  =  (Ayn/2)max 
bo  that 


co8<yAf3)max»  3in<yAf3>maxS{yAf3>max  (12> 

the  expressions  for  the  components  of  dl  and  S  in  the  X',  Y\  Z'  co¬ 
ordinate  system  become 

dix 

—  =  ynf5  -  f6  +  f5VA 
Rdx 

Aty 

—  =  [f 8  -  f4  -  yn(f7  -  f3) ]  cosMm  +  ynf3  -  f4)  -  sin(>Vm  +  ynf3  -  f4) 

Rdx 

-  {[f 8  "  f4  -  yn<f7  ~  £3>lf3  3in^m  +  ynf3  “  £4>  +  f7  cos(+m  +  Vnf3  “  £4>]Va 

+  £3<£7  '  f3>  3in^m  +  yn£3  *  £4)Ya 
d/.Z 

—  =  [f8  -  f4  "  Yn<£7  -  f3>l  sinMm  +  yn£3  '  f4>  +  cos^m  +  ^n{3  "  f4> 

Rdx 

+  {[*8  "  f4  "  yn<£7  '  £3>J£3  C0SMVn  +  yn£3  “  f4*  '  f7  sin^m  +  ynf3  “  £4>}yA 

-  £3<£7  *  f3>  cos^m  +  ynf3  -  £4^Xa 
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Sx  tr 

R 

Sy 

—  =  Yp  +  x  sin(4im  +  ynf3  -  £4)  +  xf3  cos(*\>m  +  ynf3  -  f4)yA 
R 

SZ 

- =  Zp  -  x  coa(4.m  +  ynf3  -  £4)  +  xf3  ain(+m  +  ynf3  -  f4)yA  { 13) 

R 

The  Biot-Savart  law  is  used  to  compute  the  induced  velocity  due  to 
the  vorticity  representing  the  bound  circulation.  This  law  is  stated 

^  1  r  di  X  S 

w  =  — -  |  r  — —  ( H) 

4  it  Ji  S3 

where  di  is  an  infinitesimal  length  of  an  element  of  vorticity,  S  is  the 
distance  from  the  element  of  vorticity  to  the  point  where  the  velocity, 

W,  is  desired,  and  r  is  the  strength  of  the  vortex  element.  Since  only 
the  normal  component  of  W,  i.e.  ,  Wy,  is  desired,  there  is  obtained 
from  Eq.  14 

1  f  _  SX^Z  -  Szdix 

Wy  =  —  I  T  - - -  (15) 

4  ir  Jl  S3 

The  strength,  r,  of  the  vortex  element  as  a  function  of  the  spanwise 
and  chordwise  distribution  of  bound  circulation  is  determined  as  follows. 
Specify  a  spanwise  distribution  of  bound  circulation,  r,  extending  from 
the  hub,  x  =  x^,  to  the  tip,  x  =  1,  as  illustrated  in  Fig.  5.  Assuming 
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that  the  shape  of  the  chordwise  loading  is  independent  of  3panwise  posi¬ 
tion,  it  is  necessary  to  specify  it  at  only  one  position.  Thus,  specify  a 
chordwise  distribution  of  bound  circulation  at  the  spanwise  position 
r  =  rmax,  as  shown  in  Fig.  6. 


FIG.  6. 


The  strength,  r  ,  may  then  be  written 


-  r  drmax 

F  =  - - -  dy 

^max  dY 

Nondimensionalizing  the  bound  circulation  using  the  propeller  diameter, 
D,  and  the  forward  velocity  of  the  propeller  relative  to  the  undisturbed 
fluid,  v8, 

r 

G  = - 

ttDv8 

the  expression  for  F  becomes 

F  G  dGmax 

- -  S - dy  ( 16) 

wDvs  Gmax  dY 

Restricting  tlie  chordwise  distribution  of  bound  circulation  to  one  that 
can  be  made  up  of  straight  line  segments,  as  discussed  previously, 

Eq.  16  may  be  written  for  each  spanwise  strip 


r 


ttDvs 


•'max 


dG 


max  ] 


dy 


<Y=Yn 


/d^G  max\ 

+i— L 


y=yn 


dY* 


(17) 


Equation  2  allows  the  components  of  df  and  S  in  the  X1,  Y1,  Z1  system 
to  be  determined  from  their  components  in  the  X1,  Y1,  Z1  system, 
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given  by  Eq.  13.  Substituting  these  X',  7',  Z'  components  into  Eq.  15  and 
using  Eq.  17,  the  expression  for  the  normal  component  of  induced  velocity 
due  to  the  bound  circulation  becomes 


(Wy)b  1 


v8 


-l-t  i  f  / 

2  m=  1  n=  1  •'Xh  Gmax  J- 


AVn/2 

^yn/2 


A0Hn+(A0Kn+A1Hn)yA+{A1Kn+A2Hn)y2+(A2Kn+A3Hn)y|+A3Kny^ 


(a  +  byA  +  cyJ)3/2 


(18) 


dyAdx 


where 


fdG 


Hn  = 


Kn  = 


max 


dy  / y=yn 


fd^G 


max 


dy‘ 


'y=yn 


a  =  2YpXsin(4)m+ynf3-f4)-2ZpXC08(4Jm+ynf3-f4)+(Xp-ynfi+f2)2+Yp+Zp+x2 
b  =  2Ypxf3cos(ijjm+ynf3-f4)+2Z  xf3sin(i|<m+ynf3-f4)-2f1(Xp-ynf1+f2) 


2f2  ^  r2 

c  =  x  f3  +  f  x 


■^0  =  {[f8-f4-yn(f7-f3)](Xp-ynfl+f2)sin(+m+ynf3-f4-ap)+Yp(yn£5'fG)sinaP 
+[xp-ynfl+f2+x(ynf5-f6)]c°s(‘lJm+ynf3-f4-“p)-zp(ynf5-f6)COSQp}sinPp 
-[{Yp[f8-f4-yn  f7-f3)l+zp}sin(+m+ynf3-f4)-{Zptf8-f4-yn(f7-f3)]-Yp} 

003(4^ +ynf3-f4)+x[f8-f4-yn(f7-f3)]]cospp 
A1  =  -[{f7(Xp-ynfl+f2)+fl[f8-f4-yn(f7-f3)l+xf3(ynf5-f6)}sin(+m+ynf3-f4-Qp) 
-{[f8-f4-yn(f7-£3)]f3(xp-ynfl+f2)-fl+5cf5}cos(4Jm+ynf3"f4-Qp) 

-  Ypf5  sin  ap  +  Zpf5  cos  apj  sin  (3p 
-[{Yp[f8-f4-yn(f7-f3)]£3+Zpf7}cos^m+ynf3-f4) 

+  {zp[f8-f4-yn{f7-f3)]f3-Ypf7}sin(4;m+ynf3-f4)-x(f7-f3)]cos(3p 
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A2  =  -({[f8-f4-yn{£7'f3)]flf3+(£7”f3)f3(Xp-ynfl+f2)}cos(4'rn+ynf3"f4-Qp) 

“(£1£7  -  x£3£5)sin(4)m  +  Yn£3  -  £4  -  Qp)j  s*n|3p 

+{£3(£7-£3)[YpC0S(4Jm+yn£3"£4)+Zp8in(v|jm+ynf3-£4)] 
-x[£8-£4-yn{f7-f3)]f3]co8  pp 

a3  =  £lf3(f7-£3)C0S('}'m+yn£3-£4-Qp)8inPp+xf3(£7-£3)C08  Pp 

The  integration  with  respect  to  yA  in  Eq.  18  may  be  carried  out  analytically. 
The  integration  with  respect  to  x  must  be  carried  out  numerically.  Perform¬ 
ing  the  integration  with  respect  to  yA,  and  since  the  numerical  integration 
with  respect  to  x  cannot  be  carried  through  the  singularity,  Eq.  18  becomes 


<W?>B  1  A  A  /  f*P‘*  G 


■---i  ±{f 


Qdx  + 


^max 


f 


G  \  (AWyJg 

- Qdx  + - —  (19) 


xp+«  Gmax 


v8 


where 


c  =  0  for  m  i  1  or  yn  K  yp 
Ax 

€  =  —  for  m  =  1  and  yn  =  yD 
2  v 


#For  4ac  -  b2|  >  0.0000 1 


Q  = 


A3*n 

2c 


3 

Ya 


(a  +  byA  +  cy^)1/2 


AV JZ 
-LyJz 


5b’ 


A2xn  *  A3|^n  “  xn" 


4c 


Xa _ 

+  cv?) */2 


AlKn  +  A2 


f  3b \  3b  /15b2 

H  -  K  —  -  A,  H - K  I - 

l  n  n  2c  /  3[  n2c  n\8c2 


(a  +  byA  +  cyA) 
15b2  3a 


Ayn/2 

-Ayn/2 


8^  2c 


c(4ac  -  b2) 


(2b^  -  4ac)yA  +  Zabl^n^2 
.  (a  +  byA  +  cy£)  £/2  J.^/2 

*See  footnote  on  next  page. 
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3b' 

A  i  KTn  +  A2  Hn-Kn - A3 

2c 


3b.  /15b2  3a' 


Hn - Kn  — r - 

2c  \  80^  2c I 


:3/2 


In 


(a  +  byA  +  cyA)  1/Z  +  ’'/cyA  + 


2*Vc  J 


Ayn/2 

-Ayn/2 


2a 

2a 

lOab  \ 

A0Kn  +  A  lHn  '  A2Kn - A3 

c 

«n - Kn 

c 

4c2  /. 

4ac 

byA  +  2a 

-  b2 

Ayn/2 

2A0Hn 

4. 

2cyA+b 

,(a  +  byA  +  cy|)1//>2. 

_Ayn/2 

4ac  -  b2 

(a  +  byA  +  ey^)1^2. 

Ayn/2 

-Ayn/2 


*For  |4ac  -  b2|  ^  0.00001 


A3Kn*^I 


be* 


^yn/2 

-Ayn/2 


A2Kn  +  A3  Hn  -  Kn 


8a  \ 

b  j 


1  ’  b  \  Ayn/2 

I  wi Vi  Ly„/2 


2*/a  /  2^a  3*/a\  /  3n/zi  12a^/2' 

Al^ri  +A2l^n“  ]"A3l^n  ~ 

b  \  b  c  /  \  c  be  J 


In  *v/a  +  —  yAj 

2Va  J 


Ayn/2 

-Ayn/2 


*Under  certain  circumstances  when  x=Xp,  4ac  -  b2  =  0  (exactly  or  with¬ 
in  the  accuracy  of  the  computer).  The  solution  takes  a  different  form  at 
this  point.  Due  to  machine  accuracy,  the  test  is  made  on  a  small  finite 
quantity  rather  than  on  zero. 
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\ 


I  4a  \  /4a  3a  \  I  3a  8a2  \ 

AoK„  +  Aj K-Kh—  Ua2  Hn-.-Kn--  +A3 K—  -Kn  —  J 


‘'/a+--=yA 

2va 


/  V a  <Ja  \  /  */a  a2/2  \  I  a^/2 

Ao|Hn— --Kn —  I- Aj  Hn  — -Kn— - J  +  A2  Hn - K, 


*yn/2 

-Ayn/2 

a3/'2\ 


2c 


2c 


,3/2 


*-n 


be  / 

a*/2 
be2  i 


be 

1 


b  \2 

,  I3iy‘ 


2c2  / 

Ay n/2 

-Ayn/2 


and  where  (AWy)B  is  the  contribution  to  {Wy)B  from  the  small  region 
surrounding  the  singuiarity.  The  determination  of  (A Wylfl  follows. 


Figure  7  shows  the  small  region  surrounding  the  singularity  lying 
on  the  blade  i^l  =  0  and  having  its  center  at  the  point  yn  =  yp.  The 
location  of  points  within  this  region  are  given  by  the  nondirnensional 
coordinates  y^  and 


Singularity  region 
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Since  these  coordinates  are  nondimenaionalized  with  the  chord  and  pro¬ 
peller  radius  respectively,  they  are  very  small  compared  to  one.  The 
nondimensional  coordinates  in  the  X',  Y',  Z-  system  of  points  on  the 
blade  tj»i  =  0  are,  from  Eq.  5,  7,  and  10, 

X  =  yfj  -  h 
Y  =  -x  sin(yf3  -  f4) 

Z  =  x  cos(yf3  -  f4)  (20) 

In  the  region  surrounding  the  singularity, 

x  =  xp  +  xA 

Assuming  a  linear  variation  of  the  functions  of  x  over  the  region,  the 
expressions  for  fj  _  4  are 

fl  =  flp  +  f5p*A 
h  =  f2p  +  f6pxA 


f3  8  f3P  +  (*7p  -  f3p) 


f4  8  f4p  +  (f8p  '•  f4p)  —■ 

*re.defi,ned  E*'  10‘  Substituting  these  expressions 
f?.  Eq,2°  aad  utilizing  the  trigonometric  identities  for  the  sin  and  cos 
of  the  sum.  of  angles,  there  is  obtained 

X  =  ^lp  "  f2p  +  (yf5p  ~  f6p>xA 

Y  =  -{Xp  +  xA.)sin(yf3p  -  f4p)Cos  [y(f?p  -  f3p)  -  (f8  -  f4  )J  — 

‘  (XP  *  xA)cos(yf3p  -  f4p>sin  [y(f7p  -  f3p)  -  {fg  -  f4  )]J^ 

1  XP 

Z  =  (xp  +  xA)cos(yf3p  -  f4p)cos  [y(f7p  -  f3  )  -  (fg  -  f4 )]  — 

*P 

-  (xp  +  xA)sin(yf3p  -  f4p)sin  [y(f7p  -  f3p)  -  (f8p  -  f4p)]^  (21) 
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Restricting  {xA/xp)max  3  Ax/2Xp  to  be  very  small  compared  to  one,  so 


*sin  [y(f7p  -  f3p)  -  <f8p  -  f4p)] . 


pirn  ax 


s  Hy(£7p‘f3pMf8p-f4p)l 


xp  Jmax 


cos  [y(f7p  -  f3p)  -  (f8p  “  f4p)l  —  «  1 

'  ’‘p'max 


1  +  1  — 


*p  'max 


Equation  21  becomes 


^  -  y^ip "  ^2p  *  (y^5p  ■  ^6p)xA 
Y  a  -xp  3in<yf3p  -  f4p) 

-  { sin(yf3p  -  f4p)  +  cos(yf3p  -  f4p)[y(f7p  -  f3p)  -  (f8p  -  f4p)l}xA 
Z  =  xp  cos(yf3p  -  f4p) 

+  {cos(yf3p  -  f4p)  -  sin(yf3p  -  f4p)[y(f7p  -  f3p)  -  (f8p  -  f4p)]}xA  (23) 
Introducing  again  the  expression  for  y  within  a  strip, 

y  =  yn  +  yA  =  yp  +  yA 

using  the  trigonometric  identities  for  the-  sin  and  cos  of  the  sum  of 
angles,  and  remembering  the  restriction  that  has  been  made  concerning 
(yA^max  (8ee  Eq.  12),  Eq.  23  becomes 

x  =  yPfiP  -  f2P  +  fiPyA  +  £5PyAxA  +  (ypfsp  -  f6P)xA 

Y  =  -Xp  sinQp  -  xpf3p  cosapyA 

-  f3p {cosap  -  sin Qp[(yp  +  yA)(f7p  -  f3p)  -  (fgp  -  f4p)]}yAxA 

-  (sinQp  +  coSQp[(yp  +  yA)(f7p  -  f3p)  -  (f8p  -  f4p)l}xA 
Z  =  xp  cosop  -  xpf3p  sinapyA 

-  f3P  {sinaP  +  cosap[(yp  +  yA)(f7P  -  f3P>  ■  <f8P  -  54P)]}yAxA 

+  (cos Up  -  sin Gp[(yp  +  yA)(f7p  -  f3p)  -  (f8p  -  f4p)]}xA  (24) 


*-These  approximations  imply  certain  restrictions  on  f7  and  f8  which 
are  discussed  later. 
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Using  Eq.  24  and  9*  the  expressions  for  the  components  of  d!  and  S  in  the 
X',  Y',  Z'  coordinate  system  are  easily, ahtaiaed; 

d ix  dX'  dX 

dxA  yA=Constant  dxA  yA= Constant 

dix 

—  =  V5p  ■  f6p  +  f5pYA 

dly  dY'  dY 

dxA  yA=Constant  dxA  yA=Constant  '  ~ 

dfy 

—  =  -«map  -  co.ap((yp  +  yA)(f7p  .  £3p)  .  (f8p  .  f4p)] 

-f3p{coaap  -  «inap[(yp  +  y4)(f7p  -  £3p)  -  (f8p  -  f4p)])yA 
djtZ  dZ'  dZ 

dxA  d^  yA=Constant  dx&  yA Constant 

diz 

—  =  cosap  -  Sinap[(yp  +  yA)(f7p  -  f3p)  .  (fgp  .  f4p)] 

-f3p{sinap  +  cos  ap[(yp  +  yA)(f7p  -  f3p)  -  (f8p  -  f4p)]}yA 
SX  =  Xp  ~  =  R(Xp  -  X) 

SX 

^lp^A  “  ^5pYAxA  _  ^Vp^5p  ”  f6p)xA 

SY  =  Yp  -  Y'  =  R(Yp  -  Y) 

SY 

~  =  xpf3p  cos  apyA  +  f3p{cos  ap  -  sin ap[(yp  +  yA){f?p  -  f3p)  -  (fgp  -  f4p)]}yAxA 

+  {sinap  +  cosap[(yp  +  yA)(f7p  -  f3p)  -  (f8p  -  f4p)]}xA 
Sz  =  Z'p  ~  Z1  =  R(Zp  -  Z) 

SZ 

~  =  V3p  V*  +  f3plsin °p  *  c° 3  °p[<Yp  +  Y*){%  -  *3p>  -  (*8p  '  f4p)]}yAxA 

-{cosap  -  sin ap[(yp  +  yA)(f7p  -  f3p)  -  (f8p  -  f4p)1}x&  (25) 
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Making  a  linear  approximation  in  Eq.  17  to  the  bound  circulation 
■gqn^.tfa^ratn^V»ghfcy:cesio»>v.-  *  .  *  .  .  •  •. 


Gmax  \  Gmax/^x 


Gmax 


using  Eq.  2  to  determine  the  components  of  di  and  S  in  the  X',  Y',  Z' 
coordinate  system,  discarding  third  and  higher-order  terms  in  the 
products  of  the  small  quantities  yA  and  xA  in  the  expressions 

s&tz  -  szdlx 

S2  =  +  Sy  +  s| 

and  substituting  into  Eq.  15  yields  the  expressions  for  (AWy)g/vs  . 
(AWy)n  1  /-Ayn/2  ..Ax/2 

— —  =  -  ,  <Hn  +  KnyA)(ciyA  +  c2yi) 

V„  2  J- AYn/2  J- Ax/2 


E  +  FxA 

(aiy^  +  biyAxA  +  cix&fi 


dxAdyA  (26) 


where 


Hn  = 


dGmax 


dy  'y=yn=yp 


\  ay-  / y=yn=yp 
Cl  =  -flp  sin  pp  -  Xpf3p  cos(3p 

c2  =  f3p(flp  sinPp  +  xpf3p  008  MyP(f?P  “  £3pJ  "  (f8p  -  f4p)l 


Gmax  /v= 


Gmax 


dx 
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ai  =  'lp  +  *p'3p  ^  .  .  .... 

^1 =  2riprypf5p  ■  ^6p)  *  23cp£3p[yp(f7p  -  £3p)  ~  (£8p  ■  £4p)l 

<=1  ■  (Vp'sp  -  f6p>2  +  >  +  (yp(f7p  -  £3p>  -  <£8p  -  f4p)l2 

The  integrations  in  Eq.  26  are  straightforward  but  because  they  lead  to 
lengthy  expressions,  they  are  not  presented  here.  It  suffices  that 
Eq.  26  can  be  integrated  analytically  through  the  singularity  to  obtain 
the  Cauchy  principal  value. 

VELOCITY  DUE  TO  FREE  VORTICITY  (W^f 

The  coordinates  of  a  point  on  the  helical  sheets  in  the  rectangular 
X',  Y',  Z'  coordinate  system  were  given  as  Eq.  1  and  are  repeated 
here. 


X*  =  or  tanp 

Y*  =  -r  sin(ij;m  +  a) 

Z '  =  r  cos{ipm  +  a)  ( 1) 

Assuming  that  the  propeller  is  moderately  loaded,  the  contraction  of 
the  wake  may  be  neglected.  Hence,  each  vortex  element  of  the  free 
vorticity  sheets  may  be  considered  to  lie  on  the  helical  sheet  at  a  con¬ 
stant  radius,  as  indicated  in  Fig.  8.  Letting  df  be  an  infinitesimal 


length  of  a  vortex  line  element  making  up  the  free  vorticity  sheets,  the 
components  of  di  in  the  X1,  Y1,  Z*  system  are  obtained  from  Eq.  1  as 
follows. 


dix  dX' 
da  da 


r=Con  slant 
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dix 

-  =  r  tan  p 

da 

di  y  dY' 

da  da  r=Constant 

di  y 

- =  -r  cos(i| im  +  a) 

da 

diz  dZ' 

da  da  r=Constant 

diz 

- =  -r  sin(*j/m  +  a) 

da 


The  coordinates  in  the  X',  Y\  Z'  system  of  the  point  p  on  the  blade 
i}»|  =  0  where  the  normal  component  of  velocity  is  desired  is  designated 
using  the  subscript  p,  as  in  the  previous  section. 

xp  =  Qprp  tan  Pp 

Yp  =  -rp  sinap 

Zp  =  rp  cos Qp  (27) 

The  expression  relating  ap  to  the  position  on  the  blade,  developed  in  the 
preceding  section,  is 


Qp  =  7p^3p  ‘  ^4p 


The  X',  Y1,  Z'  components  of  the  distance,  S,  from  an  element  of 
vorticity  on  one  of  the  blades  to  the  point  p  are  obtained  from  Eq.  1 
and  27. 


Sx  =  Xp  -  ar  tanp 


SY  =  Yp  +  r  sin(tm  +  Q) 


SZ  =  Zp  '  r  cos(^m  +  Q) 


Introducing  the  nondimensional  coordinates  of  Eq.  7,  and  letting 
x  tan  P  =  \,  the  components  of  di  and  S  become 


dix 

Rda 
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d£y 

•  •  ■  .  - =  -x  cos{ijim  +  a)., 

Raa 

dlz 

- =  -x  sin(^m  +  a) 

Rda 

SX 

- =  X_  -  aX 

R 

SY 

-  =  Yp  +  x  sin(4im  +  a) 

R 

SZ 

- =  Zp  -  X  cos(^m  +  a)  (28) 

R 

where 

=  ap(Mx=Xp 
Yp  =  -Xp  sinap 
“  xp  cos  Qp 

Lifting-line  solutions  for  free-running,  unskewed,  optimum  or  lightly- 
loaded  non-optimum  propellers  yield  X=  Constant,  i.e. ,  independent  of 
the  spanwise  station  x.  This  is  not  true  for  moderately  loaded,  non¬ 
optimum,  wake-adapted  propellers,  which  are  of  interest  here.  To  be 
able  to  handle  this  variation  in  X,  the  blades  are  broken  up  into  chord- 
wise  strips  bounded  by  lines  x  =  Constant,  as  illustrated  in  Fig.  9. 
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Due  to  the  great  simpli£ication  in  the  mathematical  expressions  which 
result,  the  location  of  points  along  a  strip  are  de.« cn:ihedA^,te.i:m9*M^thAv....T. 
angle  o  corresponding  to  the  center  of  the  strip.  This  causes  the  strips 
not  to  fit  the  blade  shape  exactly,  as  illustrated  in  Fig.  9  by  the  fact  that 
the  ends  of  the  strips  do  not  coincide  with  the  leading  edge.  However, 
by  using  a  fairly  large  number  of  strips,  the  blade  shape  can  be  repre¬ 
sented  quite  accurately.  For  each  strip  there  may  be  written 

x  =  xn  +  xA  (29> 

The  subscript  n  takes  on  values  of  1,  2,  3,  , . . ,  i,  where  i  is  the  num¬ 
ber  of  strips.  With  the  blade  broken  up  into  a  fairly  large  number  of 
strips,  the  variation  in  X  can  be  represented  with  sufficient  accuracy 
by  a  linear  approximation  across  each  strip. 


Substituting  Eq.  29  and  30  in  Eq.  28,  the  expressions  for  the  components 
of  di  and  S  become 


d<X 

Rda 


MM 


x=xn 


+ 


d|y 

- -  -xn  cos(4>m  +  a)  -  cos(4im  +  a)xA 

Rda 


dljr 

- =  ^  sin(4>m  +  a)  -  sin(4>m  +  a)xA 

Rda 


Sx  /  d\  \ 

- =  Xp  -  a(Mx=xn  "  Q  ( — I  XA 

R  \dx/x=xn 

SY 

- =  Yo  +  xn  sin^m  +  Q)  +  sin(4»m  +  q)xa 

R 


SZ 

-  =  Zp  “  xn  cos^m  +  a)  -  cos<tm  +  <*)xA  (31) 

R 


The  Biot-Savart  law,  Eq.  14,  is  used  to  compute  the  induced  ve¬ 
locities  from  the  free  vorticity.  Since  only  the  normal  component  of 
velocity  is  desired,  Eq.  15  gives  the  desired  relationship.  The  strength, 
F,  of  the  vortex  element  as  a  function  of  the  spanwise  and  chordwise 
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distribution  of  free  vorticity  is  determined  as  follows.  The  value  of 
I  at  thc-'trarKirg  edgv?'oHsttr''»ratfe,"  P'T^fFfgf &Y$  itr- terms- of  the  bound 
circulation  T  is  given  by, 

dr 

rTjr  =  - —  dx 
TE  d* 

A  factor  k  is  defined  as  the  ratio  of  the  strength  P  at  any  chordwise 
station  to  the  strength  at  the  trailing  edge,  FTE. 

r 

k  =_ - 

r  TE 

Introducing  once  again  the  nondimensional  bound  circulation  G  =  r/irDvs 
the  expression  for  T  is 


r  dG 

- =  -k — dx  (32) 

ttDv3  dx 

Considering  an  infinitesimal  area  on  the  blade  bounded  by  lines 
y  =  Constant  and  x  =  Constant  and  utilizing  the  fact  that  the  sum  of  the 
bound  circulation  and  free  vorticity  entering  the  area  must  equal  the 
sum  of  the  bound  circulation  and  free  vorticity  leaving  the  area,  the 
expression,  illustrated  in  Fig.  10,  for  k  may  be  determined. 


FIG.  10. 
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1  /*y  dGma3C. 

For  0  <  y  ^  1  k  = -  |  - dy  For  y  >  1  k  =  1 

Gmax  Jo  dy 

This  simple  relationship  between1  k  and  y,  which  is  independent  of  the 
spanwise  station,  results  from  having  chosen  the  elements  of  bound 
circulation  along;  lines  y  =  Constant.  To  obtain  k  as  a  function  of  a,  the 
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c  =  1  -  q2n£ 

Ao  =  -{(Xp  -  aMn)xn  sin(4im  +  a  -  ap) 

"  MJxn  C03(+m  +  Q  -  Qp>  +  Yp  sinap  -  Zp  c°8Qp]}9inPp 

+  xjYp  8in^m  +  a)  '  zp  coa(4'm  +  a)  +  xjcos  pp 

A1  =  ’K<Xp  ‘  aMn>  "  sinW*m  +  Q  “  Qp) 

-  (Mn  +  Nnxn)cos(4)m  +  a  -  ap)  -Nn(Yp  sinap  -  Zp  co8ap)}sinpp 

+  [Yp  sin(4*m  +  a)  -  Zp  ccs(i]/m  +  a)  +  2xn]cos  pp 
A2  =  Nn[a  sin(4>m  +  a  -  ap)  +  cos(4<m  +  a  -  Qp)]sinpp  +  cospp 

and  where 


=  (^x=xn 

I  d\\ 


Nn  = 


idx 


X=Xr 


The  integration  with  respect  to  xA  in  Eq,  34  may  be  carried  out  analyti¬ 
cally.  The  integration  with  respect  to  a  must  be  carried  out  numeri¬ 
cally.  Performing  the  integration  with  respect  to  xA,  and  since  the 
numerical  integration  with  respect  to  a  cannot  be  carried  through  the 
singularity,  Eq.  34  becomes 


(W  y)f 

vs 

where 


,  g  i 

1  V'  V'  /  faP'(  /*au 

=  -—  2^  [  I  kQda  +  |  kQda 

2  m=l  n=l  \  JaLE  Ja p+t 


(AWy)p 

+ - (35) 


e  =  0  for  m  /  1  or  *n  /  xp 


ha 

e  =  —  for  m  =  1  and  xn  =  Xp 
2 
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*Under  certain  circumstances  when  a  =  Qp,  4ac-b^  =  0.  The  solu¬ 
tion  takes  a  different  form  at  this  point.  Due  to  machine  accuracy,  the 
test  is  made  on  a  small  finite  quantity  rather  than  on  zero. 
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4a 

4a 

3a\ 

A0Jn  +  A1 

fn  *  Jn  - 

-  a2 

rn-T 

-  Jn- 

b 

b 

c  / 

c 


\/a  +  • 


2'>/a 


Axn/2 


-Axn/2 


/  *'/a  Va  \ 

M'-T-1"*  h'1' 


*ya"  a^/  2 


-J. 


\ 


t3/2 


iV2\l 


Mn2c  “n  be  / 


+  A2  rn“ - J 


be 


"  2c2 


j 


AXn/2 


-Axn/2 


and  where  (AWy)F  23  the  contribution  to  (Wy)p  from  the  small  region 
surrounding  the  singularity.  The  determination  of  (AW y)p  will  follow 
shortly.  The  value  of  a  at  the  leading  edge,  a^g,  follows  from  the 
work  in  the  previous  section. 


aLE  =  -tf4)x=xn 

The  value  of  the  upper  limit  of  integration,  au,  is  theoretically  in¬ 
finity;  however,  as  is  discussed  subsequently,  the  integration  maybe 
truncated  at  a  rather  small  value  of  qu.  Assuming  that  the  integration 
is  carried  out  q  propeller  diameters  (measured  axially)  behind  the 
point  p  where  the  velocity  is  desired,  the  expression  for  qu  is 

2cl  +  ap(^X=Xp 

Q  u  = - E 

Wx=xn 

This  follows  from  the  first  of  Eq.  1. 

Since  the  lifting-surface  solution  starts  from  the  lifting-line  solu¬ 
tion,  the  angularity  of  the  flow  due  to  the  lifting-line  solution  has 
already  been  accounted  for.  Hence,  only  the  difference  between  the 
normal  component  of  velocity  on  the  lifting  surface  and  the  normal 
component  of  velocity  at  the  lifting  line  from  the  lifting-line  solution 
is  desired.  It  is  of  course  understood  that  both  these  velocities  are 
determined  at  the  same  radius.  For  free-running,  optimum  or  lightly 
loaded  non-optimum  propellers  where  the  helical  sheets  form  a  true 
helix,  i.e.,  X.  =  x  tan  p  =  Constant,  it  i s  necessary  in  determining  this 
difference  to  choose  a  value  of  q  which  assures  only  that  the  integrations 
are  carried  cut  to  a  point  behind  the  trailing  edge.  This  is  true  because 
for  the  true  helix  the  integration  from  q  diameters  behind  any  point  (at 
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a  given  radius)  to  infinity  will  be  the  same,  so  that  upon  taking  the  dif¬ 
ference  thi a. contribution  cancels  out  exactly. _ This  fact  hasubeen  taken 

advantage  of  by  several  investigators  in  calculating  lifting- surface 
corrections  for  free- running,  optimum,  or  lightly  loaded  non-optimum 
propellers.  For  all  other  propellers,  such  as  the  moderately  loaded, 
non-optimum,  wake -adapted  propellers  of  interest  here,  the  helical 
sheets  do  not  form  a  true  helix  and  this  canceling  effect  does  not  occur. 
Luckily  for  these  cases,  the  difference  between  the  lifting-surface 
solution  and  the  lifting-line  solution  is  strongly  dependent  on  q  only  for 
small  values  of  q.  Thus  it  is  necessary  to  carry  out  the  integration 
only  a  few  propeller  diameters  behind  the  two  points  when  determining 
this  difference.  A  great  saving  in  computer  time  results  from  this 
approach,  since  numerical  integrations  to  great  distances  behind  the 
propeller  are  not  needed. 

The  normal  component  of  velocity  at  the  lifting  line  from  a  lifting¬ 
line  solution  where  the  integration  is  carried  out  q  propeller  diameters 
behind  the  lifting  line  is  needed  to  determine  the  difference  discussed 
above.  This  is  determined  from  Eq.  35  by  putting  ap  =  0,  k  =  1,  and 
ignoring  the  integral  from  ole  to  cip  -  c.  The  expression  for  (AWY)F 
takes  a  different  form  for  the  lifting-line  solution,  as  will  be  indicated 
later. 

Figure  11  shows  the  small  region  surrounding  the  singularity.  This 
region  lies  on  the  blade  =  0  and  has  its  center  at  the  point  xn  =  Xp. 


The  location  of  points  within  this  region  is  given  by  the  nondimension^l 
coordinate  x^  and  the  angle  q^.  Since  x^  is  nondimensionalized  with 
respect  to  the  propeller  radius,  it  is  very  small  compared  to  one.  Also, 
is  restricted  to  be  very  small  compared  to  one.  The  nondimensional 
coordinates  in  the  X1,  Y',  Z'  system  of  points  on  the  blade  ijjj  =  0  are, 
from  Eq.  1  and  7 
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X  =  ax  tan  (3  =  aX. 

Y  =  -x  sin  a 

Z  =  ::  cos  a  (36) 

In  the  region  surrounding  the  singularity, 


x  =  xn  +  =  xp  +  xA ,  a  =  ap  +  aA 

Substituting  these  expressions  and  Eq.  30  into  Eq.  36,  and  making  use  of 
the  trigonometric  identities  for  the  sin  and  cos  of  the  sum  of  angles, 


X  *  (ap  +  aA) 


Id \\ 

(Mx=Xp  +  — ) 


dx/x=xp 


*A 


Y  =  -{xp  +  xA)(sinap  cosaA  +  cos  ap  sinaA) 
Z  =  (Xp  +  xA)(coscip  cosnA  -  sinQp  sinaA) 


(37) 


Since  (aA)max  =  Aa/2  is  restricted  to  be  very  small  compared  to  one 
(in  radian  measure), 


cosaA  S  1,  sinaA  S  aA 

Putting  these  approximations  in  Eq.  37,  there  is  obtained 


(d\ ' 


fd\) 


X  =  Qp(X)x=XP  +  Qp\dx/x=xp  **  +  \dx)x=xp  XA°A  +  (MX=XP  °A 

Y  =  -xp  sinap  -  sinapxA  -  cosQpxAQA  -  xp  cosapQA 

Z  =  Xp  cos  ap  +  cosapXA  -  sindpXAnA  -  Xp  sindpdA  (38) 

Using  Eq.  38  and  the  expressions  on  page  24  for  XpYpZp,  the  expressions 
for  the  components  of  df  and  S  in  the  X1,  Y1,  Z‘  coordinate  system  are 
easily  obtained. 


dfx 

ddA 


dX' 

ddA 


dX 

=  R - 

x^aConstant  daA 


x^=Constant 


dix 

RdoA 


+ 
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d!  Y  dY' 


daA 


daA 


xA=Constant 


dY 

R - 

da^ 


xA=Constant 


diy 
RdaA 

dlz  dZ1 
daA  daA 

diz 


-xp  cosap  -  cosapxA 


dZ 


=  R 


x*=Constant 


xA=Constant 


=  -xp  sin ap  -  smapxA 


RdaA 

sx  =  Xj,  -  X'  =  R(Xp  -  X) 


Sx 

/ 

d\\ 

[d\ 

-  = 

"qd 

— 

i  XA 

-  [  — 

R 

P\ 

dx , 

X=Xp 

\  dx 

sY  = 

V 

Y' 

=  R(Yp  • 

■  Y) 

Sy 

R 

sin  a 

pxA 

+  cosap 

xAqA 

SZ  = 

ZP- 

Z' 

=  R(zp  ■ 

•  7 A 

=  -COSQpxA  +  sinQpxAQA  +  Xp  sinapaA 


Making  a  linear  approximation  in  Eq.  33  to  k  across  the  singularity 
region, 


k  =  <k>a=ap  + 


dk\ 

da  / a=ar 


using  Eq,  2  to  determine  the  components  of  df  and  S  in  the  X',  Y\  Z* 
coordinate  system,  discarding  third  and  higher-order  terms  in  the 
products  of  the  small  quantities  xa  and  ua  in  the  expression 
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S“  =  Sjr-  +  Sy  +  % 


and  substituting  into  Eq.  15  yields  the  expression  for  (iWy)p. 

*Aa/2 
-Aa/2 

k0  +  kl  QA 


(AWy)p  1  n/2  ? 

dn  +  Jn^HC^  +  C2x£) 

2  *'-Axn/2 


£ 


where 


In  = 


Jn  = 


+  blxAQA  +  c1qA>3/2 


f  dG  \ 
idx  /x=Xn= 


daAdxA 


xn=xp 


'd2Gi 


dx2  /  x=xn=xp 


k0  =  (k)a=a 


k  !  =  I  — 


P 

dk\ 

da  /q=qt 


and  where 


Ci  =  Mn  sin  Pp  +  xp  cosPp 

Cz  =  Nn  sin  pp  +  cospp 

.  2..  2 

al  =  1  "  apNn 


bl  =  2*pMnNn 
ci  =  M,y  +  x^ 


=  d)x=xn=Xp 
dx/x=xn=xp 


(39) 
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For  the  lifting-line  solution,  where  the  point  p  is  at  the  front  edge  of  the 
helicaLshsefcan(Lther.e„i*.na  vaxiatiacuiirthA  strength,  of  the^-shAeLwith.  - 
respect  to  a,  the  lower  limit  of  integration  with  respect  to  aA  should  be 
zero  instead  of  -Aa/2,  and  the  values  of  ko  and  kj  should  be  one  and 
zero,  respectively.  The  integrations  in  Eq.  39  can  be  carried  out  ana¬ 
lytically  through  the  singularity  to  obtain  the  Cauchy  principal  value. 
These  integrations  lead  to  lengthy  expressions  that  are  not  presented 
here. 

VELOCITY  DUE  TO  BLADE  THICKNESS  (Wy)T 

In  calculating  the  effect  of  blade  thickness,  the  blades  are  broken 
up  into  spanwise  strips  like  those  used  in  the  section  on  bound  circula¬ 
tion.  Hence,  many  expressions  take  the  same  form  and  are  taken 
directly  from  that  section. 

The  components  of  the  distance,  S,  from  a  point  on  one  of  the  blades 
to  the  point  p  where  the  velocity  is  desired  are  given  by  Eq.  13  and  are 
repeated  here. 


- =  XP  -  ynfi  +  f2  -  fiy*  • 

SY 

- =  Y  +  x  sinM>m  +  ynf3  -  f4,  +  xf3  cos(«|»  +  ynf3  -  f4)yA 

R 

SZ 

-  =  Z  -  x  cos(4;  +  ynf3  -  f4)  +  xf3  sin(^  +  y  fj  -  f4)yA  ( 13) 

R  ^ 

Since  the  velocity  at  point  p  caused  by  a  source  (representing  the  blade 
thickness)  on  one  of  die  blades  is  desired,  the  source  potential  may  be 
written 


q  1  q  1 

4tt  S  4tt  ^S^  r  S^  + 
where  q  is  the.  source  strength. 

The  components  of  the  induced  velocity  at  point  p  in  the  X1,  Y1,  Z' 
directions  from  the  source  are  given  by 

q  Sx 

asx  4tt  (sx  +  s|  +  s|)3/2 

94>  q  Sy 

3Sy  '  Y  4  it  (Sx  +  Sy  +  s|)3/2 
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8<t>  n  sz 

3SZ  z  4ir  (sj  +  s^.  +  s|)V2 

Equation  2  allows  the  normal  component  of  velocity,  Wy,  to  be  de¬ 
termined  from  the  preceding  result.  The  normal  component  of  velocity 
at  point  p  due  to  a  source  of  infinitesimal  strength  on  one  of  the  blades' 
is  then  written 


dq  Sjf  cosPp  +  Sy  cos  Op  sin(3p  +  sinap  sinPp 
4ir  (S^  +  S  y  +  s|)3/2 


dWy  =  — 


(40) 


The  infinitesimal  source  strength  can  be  written  as  a  product  of  a  sur¬ 
face  source  density  on  the  helical  sheet,  p,  and  an  infinitesimal  area, 
dA,  of  the  sheet. 


dr)  =  pdA  (41) 

The  3ource  density  can  be  related  to  the  blade-section  thickness  distri¬ 
bution  as  is  done  in  slender-airfoil  theory  and  illustrated  in  Fig.  12. 

dz  * 

p  s  2V  —  (42) 

dy 


Normal  velocity  at  sheet 
V 


Sheet  of  sources 
of  local  density, p 


=  —  -1. 
C  C 


.  dz 

pS2V  — 
dy 


FIG.  12. 


*Since  dz/dy  goes  to  infinity  as  y  -*■  0  for  essentially  all  thickness 
forms,  the  nose  shape  is  modified  slightly  to  retain  a  finite  dz/dy  at 
y  =  0. 
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The  relative  velocity  V  of  the  blade  section  and  the  fluid  is  given  approxi¬ 
mately  for  a  propeller  by 

V  =  \j(v  +  Wa)2  +  (ur  -  Wt)2 


where  v  is  the  local  axial  inflow  velocity  to  the  propeller,  w  is  the  angu¬ 
lar  velocity  of  the  propeller,  and  Wa  and  are  the  axial  and  tangential 
components  of  induced  velocity  from  the  lifting-line  solution.  Assuming 
that  the  thickness  distribution  function,  z  =  f(y),  is  the  same  at  all  span- 
wise  stations  except  for  a  multiplier  that  is  a  function  of  x  only,!  it  is 
necessary  to  specify  it  at  only  one  station.  Hence, 


dz 

/dz\ 

—  = 

I— 1 

dy 

\dy/Ref 

(43) 


where  r^  is  the  ratio  of  the  percent  thickness  of  the  blade  section  at  any 
spanwise  station  to  the  percent  thickness  at  the  reference  station.  The 
expression  for  the  infinitesimal  area,  dA,  of  the  helical  sheets  in  terms 
of  the  nondimensional  coordinates  y  and  x  can  be  determined  by  the  use 
of  Eq.  4,  5,  7,  and  10  and  relationships  involving  direction  cosines.  The 
expression  is 

dA  =  RC\Z{(yf5  -  ffc)cos  (3  +  [f8  -  f4  -  y(f7  -  f3)]sinp},i  +  1  dydx 

Since  the  variation  in  y  across  a  strip  is  not  large,  a  fairly  accurate 
approximate  expression  for  dA  may  be  obtained  by  setting  y  equal  to 
the  value  at  the  center  of  the  strip,  yn. 

dA  =  RCFndydx  (44) 

where  _ _ 

Fn  =  V{< Yrvf 5  “  f6>cosP  +  f8  "  f4  '  Yn([7  '  f3>lsinP}2  +  1 

Substituting  Eq.4l,  42,  43,  and  44  into  Eq.40,  the  normal  component  of 
velocity  at  point  p  due  to  the  thickness  effect  of  an  infinitesimal  area  on 
one  of  the  blades  becomes 


*This  type  of  relation  holds  approximately  for  any  of  the  NACA 
basic  thickness-form  groups.  For  example,  the  z  =  f(Y)  for  the  NACA 
65-008  form  can  be  obtained  approximately  from  that  for  the  NACA  65- 
010  by  multiplying  by  8/l0, 
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a  =  2Ypx  sin{ij/m  +  yn£3  -  f4)  -  2Zpx  cos(4»m  +  ynf3  -  £4) 

+  <XP  -  Vi  +  (2)2  +  +  zp +  *2 

b  =  2Ypxf3  +  ^3  -  f4> 

+  2Zpxf3  sin(i}jm  +  ynf3  -  £4)  -  2£j(Xp  -  y,^  +  f2) 

c  =  x2f2  +  f  ] 

The  integration  with  respect  to  yA  in  Eq.  47  can  be  carried  out  analyti¬ 
cally.  The  integration  with  respect  to  x  must  be  carried  out  numeri¬ 
cally.  Carrying  out  the  integration  with  respect  to  yA  and  since  the 
numerical  integration  with  respect  to  x  cannot  be  carried  through  the 
singularity,  Eq.47  becomes 


(Wy)^i  1 


•- i  i ( r* w—r 

*  "l  ISm  '  ■/xP+< 


where 


«  =  0  for  m  /  1  or  yn  /  yp 
Ax 

«  = -  for  m  =  1  and  yn  =  yp 

2 


1  \  (AWy)'j1 

foFnQdx  + - 


*For  |4ac  -  b2|  >  0.00001 
Q  = 


2A0Sn 

2cyA  +  b 

4ac  -  b2 

.(a  +  byA  +  cy2)  l/2 

2(A0Tn  +  AlSn) 


4ac  -  b2 


AlTn 

c(4ac  -  b2) 


*y  n/2 

-Ayn/2 

byA  +  2a  lAyn/2 
(a  +  byA  +  cy2) 1/2 


(2b2  -  4ac)y^  +  2ab 


(a  +  by  +  cy2)  l/z 


-Ayn/2 

AYn/2 

-AYn/2 


AlTn  [ 


In 


(a  +  byA  +  cy^)  ^/z  +  'Tcy  + 


Z*fc 


Ayn/2 

-Ayn/2 


*Under  certain  circumstances  when  x  =  xp,  4ac  -  b2  =  0  (exactly  or 
within  the  accuracy  of  the  computer).  The  solution  takes  a  different 
form  at  this  point.  Due  to  machine  accuracy,  the  test  is  made  on  a 
small  finite  quantity  rather  than  on  zero. 
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using  the  value  of  Fn  at  the  center  of  the  region, 

Fn  =  (Fn)x=xp  =  FP 

Yn=yP 

discarding  third  and  higher-order  terms  in  the  products  of  the  small 
quantities  yA  and  xA  in  the  expressions 

(dz\ 

—  (sx  cos  Pp  +  Sy  COS  Qp  sin  Pp  +  Sz  sin  a  sin  Pp) 
\dy/Ref 

s2  =  +  s| 

substituting  into  Eq.45  and  indicating  the  integration  over  the  region, 
the  expression  for  (4Wy)^  is 


(AW  y)y 
vs 


1  /*Ax/2 

-  F  I  (A  +  Bxa) 

w  •'-Ax/2 


coxa  +  cixaYa 

- - - T-rjr  dyAdxA 

(alXA  +  Va^A  +  C1Ya)  7 


(48) 


where 

Fp  =  V{(Ypf5p  -  f6p)  cos  Pp  +  [f8p  -  f4p  -  yp(f7P  -  f3p)lsinPp}2  +  1 

A  =  ^O^XSXp 

B  (df°l 

\  dx  /  x=Xp 

C0  =  Sn(typ(f7p  “  f3p)  '  (f8p  -  f4p)lsinPp  *  <ypf5p  “  f6p>COS(3p} 

^  1  =  ^n{typ^7p  "  "  ^8p  “  ^4p)ls'nPp  "  (Ypf5p  "  ^6p)  cos  Pp} 

*  Sn^7p  s*nPp  "  %p  cosPp) 

SnsH 

\d>  'y=yn=yp 


'd2zi 


Tn  = 


idy2/y=yn=yp 
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ai  8  (yPf5P  ■  f6P)2  +  1  +  CyP(f7P  -  f3P)  -  (f8p  •  f4p)]2 
bj  ~  ^lp(yp*5p  “  *6pf  ^pt3piyp(i7p  "  ^3pJ  “  I r8p  "  ^4pl l 


The  integration  in  Eq, 48  can  be  carried  out  analytically  through  the 
singularity  to  determine  the  Cauchy  principal  value.  Due  to  the  lengthy 
expression  that  results,  this  operation  is  not  shown  here. 


DETERMINATION  OF  CAMBER  LINES 
FROM  NORMAL  COMPONENT  OF  INDUCED  VELOCITY 

The  total  normal  component  of  induced  velocity  desired  for  comput¬ 
ing  the  camber  lines  is 

W y  =  (Wy)B  +  (WY>F  "  (WY)FLL  +  (wy)T 

where  (Wy^FLL  *s  the  normal  component  of  induced  velocity  at  the 
lifting  line  from  the  lifting-line  solution.  This  component  must  be  sub¬ 
tracted  since  the  angularity  of  the  flow  due  to  the  lifting- line  solution 
has  already  been  accounted  for.  This  point  was  discussed  in  detail  in 
the  section  on  free  vorticity.  The  nomenclature  used  in  describing  the 
camber  lines  is  given  in  Fig.  13. 


FIG.  13. 

The  camber  lines  are  approximated  by  a  power  series  of  j  4-  1  terms, 
h  =  ag  +  ajy  +  a£y2  +  a3y3  +  •  •  *  +  ajyj 
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The  slope  of  the  camber  line  is  related  to  the  normal  component  of  in¬ 
duced-  velocityoy  ' 

dh  W  y 

dy  V 

Hence, 

W  Y  2 

- =  ai  +  2a  >y  +  3a?v  +  •  •  •  +  ja;yJ 

V  J 

If  Wy  is  determined  at  j  values  of  y  across  the  blade,  there  result 
j  linear  equations  in  the  j  unknowns,  aj.  Upon  solving  for  the  aj,  the 
camber  offset,  Ah,  and  the  variation  of  the  pitch  angle  from  the  lifting 
line  value,  y,  are  determined  from 


i=  1 


Ah  =  (y  -  l)a  +  (y  -  l)aj  +  (y2  -  l)a2  +  •  •  •  +  (y*  -  l)aj  (49) 


COMPUTER  PROGRAMS 

The  foregoing  equations  for  the  solution  of  the  induced  velocities  and 
the  camber  lines  were  programmed  on  an  IBM  7090  computer.  To  fa¬ 
cilitate  program  checkout,  four  separate  programs  were  written:  (1) 
solution  of  induced  velocities  due  to  bound  circulation,  (2)  solution  of 
induced  velocities  due  to  free  vorticity,  (3)  solution  of  induced  velocities 
due  to  blade  thickness,  and  (4)  solution  of  camber  lines  and  angles  of 
attack  from  induced  velocities.  The  induction  factor  lifting-line  solution 
of  Lerbs  (Ref.  2)  was  also  programmed  on  the  7090  computer  because  it 
serves  as  the  starting  point  of  the  lifting- surface  solution.  These  five 
programs  were  thoroughly  checked  out  against  hand  calculations  and 
special  cases  for  which  exact  solutions  can  be  obtained.  Some  of  the 
results  and  some  typical  design  calculation-  for  a  wake-adapted  pro¬ 
peller  will  be  given  in  a  future  report  on  numerical  results  from  the 
lifting-surface  theorv. 

Simpson's  rule  was  used  in  these  programs  when  numerical  integra¬ 
tions  were  needed.  Although  there  are  methods  superior  to  Simpson's 
rule,  it  appeared  more  advantageous  to  gain  accuracy  by  using  more 
points  in  the  numerical  integration  than  to  use  a  more  complicated  inte¬ 
gration  technique. 

Examination  of  the  foregoing  mathematical  developments  shows  that 
first  and  second  derivatives  of  various  functions  arc  needed  as  inputs 
to  the  solution.  Normally,  these  functions  are  known  only  as  a  table  of 
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discrete  values  and  determination  of  meaningful  third  or  higher  deriva¬ 
tives  is  very  difficult.  Consequently,  the  approach  taken. * 
ceding  mathematical  developments  was  adopted,  e.  g.  ,  low-order  ap¬ 
proximations  (involving  at  most  second  derivatives)  were  used  across 
the  strips  and  accuracy  obtained  by  using  many  strips.  The  needed 
first  and  second  derivatives  could  then  be  determined  through  a  7090 
computer  program  existing  at  NOTS  that  allows  a  sliding  polynomial 
fit  by  least  squares  to  the  discrete  input  data  and  outputs  first  and 
second  derivatives. 


RESTRICTION  ON  BLADE  SHAPE 

Blade  shapes  having  rounded  tips,  as  illustrated  in  Fig.  14,  cannot 
be  handled  by  this  method  because  of  the  behavior  of  certain  functions 
that  appear  in  the  numerical  integrations. 

Allowable  blade  tip  shape 

o 

\  /  Round  tip  shape 

\  /  which  cannot  be 

\  /  handled  by  design 

I  (  method 

FIG.  14.  - 


d[2(C/D)  sin  p] 

,  - - - - -  -  00 

dx 

d[2(L/D)  sin  p] 

- - - -  -  00 

dx 

d[2(C/D)  cos  p] 

- -  CO 

dx 

d[2(.L/D)  cos  p] 

- — -  oo 

dx 


For  rounded  tips, 

*5  = 

f6  = 

h  = 

f8  = 
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Also,  the  approximations  of  Eq.  22  can  be  made  valid  only  if  (7  and  fs 
do  oat.be, come  excessively  large.  Hence,  this  method  is  restricted  to 
blade  shapes  that  have  a  finite  chord  at  the  tip,  as  shown  in  Fig.  14. 
Since  such  a  modification  to  the  more  usual  rounded  tip  can  result  in 
at  most  a  very  small  loss  in  propeller  efficiency,  no  significant  loss 
in  the  method's  usefulness  is  caused  by  this  blade-shape  restriction. 


DISCUSSION 

The  approach  taken  in  the  preceding  work  has  been  to  obtain  a  satis¬ 
factory  engineering  solution  to  a  difficult  problem  rather  than  to  attempt 
an  elegant  development  from  a  mathematical  point  of  view.  Therefore, 
this  method  must  be  used  with  a  certain  amount  of  engineering  judgment 
in  such  things  as  the  number  and  width  of  the  strips  into  which  the  blades 
are  divided  and  the  size  of  the  singularity  region.  The  check  solutions 
and  design  calculations  that  have  been  run  to  date  provided  sufficient 
knowledge  for  these  judgments  to  be  made  now  on  a  rational  basis.  A 
further  discussion  of  this  point,  along  with  typical  values,  will  be  given 
in  a  future  report  presenting  numerical  results  from  the  lifting-surface 
theory. 

Verification  of  an  interesting  point  has  occurred  in  the  calculations 
that  have  been  run.  For  unskewed,  free- running,  optimum,  or  lightly 
loaded  non-optimum  (X  =  x  tan  p  =  Constant)  propellers  having  symmet¬ 
rical  blade  shape,  symmetrical  chordwise  loading,  and  negligible  thick¬ 
ness,  the  pitch  given  by  the  lifting- surface  solution  is  identical  to  that 
from  the  lif*ing-line  solution  (y  in  Eq.  49  equals  zero  for  all  spanwise 
stations).  This  point  has  been  referred  to  by  several  investigators. 
Hence,  the  pitch  correction  for  such  propellers  arises  solely  from 
blade  thickness. 


CONCLUSIONS  AND  FUTURE  WORK 

The  lifting- surface  design  method,  with  the  aid  of  a  high-speed 
computer,  allows  single-rotating,  wake-adapted  propellers  of  nearly 
arbitrary  shape  and  loading  to  be  designed  without  the  need  for  the 
many  assumptions  and  approximations  involved  in  the  methods  that  can 
be  carried  out  with  a  hand  calculator.  Since  the  mathematical  model 
accounts  for  nearly  all  parts  of  the  physical  propeller  system  in  a 
fairly  exact  manner,  there  is  reason  to  believe  that  propellers  designed 
by  this  lifting-surface  method  will  perform  to  specification.  A  definite 
conclusion  cannot  be  made,  of  course,  until  propellers  are  designed  and 
tested. 

In  torpedo  design,  counterrotating  propellers  are  usually  employed 
to  obtain  a  torque  balance.  In  order  to  apply  the  lifting-surface  method, 
a  lifting-line  solution  is  under  development  for  counterrotating  propellers 
that  will  serve  as  the  starting  point  for  application  of  the  lifting- surface 
solution  '•o  each  propeller.  The  effect  of  each  propeller  on  itself  is 
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determined  by  a  slightly  modified  version  of  Lerbs1  induction-factor 
method  that  allows  the  effect  of  the  interference  velocities  (those  in¬ 
duced  at  one  propeller  by  the  other)  to  be  accounted  for  in  the  relaulvc 
flow  at  the  lifting  line  and  also  allows  application  to  propellers  having 
a  finite  circulation  at  the  hub. 

The  mean  axial-induced  velocities  at  one  propeller  induced  by  the 
other  are  obtained  by  replacing  the  finitely  bladed  propeller  with  an 
infinitely  bladed  propeller  having  the  same  radial  thrust  distribution. 
The  mean  tangential- induced  velocities  at  the  rear  propeller  caused  by 
the  front  propeller  are  determined  from  Stokes'  theorem.  Continuity 
effects  due  to  hub  taper  are  accounted  for  in  an  approximate  manner. 

The  one  possibly  important  factor  in  propeller  design,  which  to 
date  has  not  been  thoroughly  investigated,  is  the  effect  of  the  hub 
boundary  condition.  Work  has  been  done  on  optimum  propellers  where 
the  hub  boundary  condition  was  satisfied  in  the  ultimate  wake;  however, 
it  has  not  been  shown  that  this  bears  any  relation  to  satisfying  the 
boundary  condition  in  the  vicinity  of  the  propeller,,  The  possibility  of 
satisfying  the  hub  boundary  condition  in  the  vicinity  of  the  propeller  by 
using  a  surface  source  density  on  the  hub  is  being  considered,  A 
method  similar  to  that  used  by  Smith  and  Hess  (Ref,  12)  would  be  used. 
Until  the  hub  boundary  condition  is  examined  in  some  such  fairly  rigor¬ 
ous  fashion,  there  will  always  remain  a  doubt  as  to  the  adequacy  of 
even  a  lifting-surface  design  method. 


Negative  Numbers  or  Illustrations:  None 
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